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THE FREE PATH IN A HIGH VELOCITY RANDOM FLIGHT PROCESS 
ASSOCIATED TO A LORENTZ GAS IN AN EXTERNAL FIELD 

ALEXANDRU HENING, DOUGLAS RIZZOLO, AND ERIC S. WAYMAN 


Abstract. We investigate the asymptotic behavior of the free path of a variable density 
random flight model in an external field as the initial velocity of the particle goes to infinity. 
The random flight models we study arise naturally as the Boltzmann-Grad limit of a random 
Lorentz gas in the presence of an external field. By analyzing the time duration of the free 
path, we obtain exact forms for the asymptotic mean and variance of the free path in terms 
of the external field and the density of scatterers. As a consequence, we obtain a diffusion 
approximation for the joint process of the particle observed at reflection times and the amount 
of time spent in free flight. 
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1. Introduction 

We study the behavior of a random flight process that arises as the Boltzmann-Grad limit of 
a random scatterer Lorentz gas in an external field. The Lorentz gas model was introduced in 
1905 by H.A. Lorentz [Lor05] as a model for the motion of electrons in metallic bodies. Since 
its introduction, the model has been widely studied by both mathematicians and physicists 
(see |Detl4| for a recent survey). In this model, a point particle travels in an array of fixed 
convex scatterers. When the particle comes in contact with a scatterer it reflects specularly. 
This variant is referred to as the “hard core” version. There is also a “soft core” version in 
which the scatterers deflect the moving particle via a potential instead of a hard boundary. Of 
the many variations on the hard core model, our study is motivated by the version in which 
an array of spherical scatterers is chosen randomly and the flight of the moving particle is 
determined by the action of an external field. The random scatterer Lorentz gas is difficult to 
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study directly even in the absence of an external field. One approach that has proven tractable, 
however, is the Boltzmann-Grad limit. The Boltzmann-Grad limit is a low density limit in 
which the number of scatterers in a fixed box goes to infinity while, at the same time, the 
size of each scatterer goes to zero in such a way that the total volume of the scatterers in the 
box goes to zero. If the centers of scatterers are placed according to a Poisson process and 
the rates are chosen appropriately, the asymptotic behavior of the moving particle is described 
by a Markovian random flight process |Gal781 |Spo78| |Spo88| . The Markovian nature of the 
Boltzmann-Grad limit is due to the following two observations: re-collisions with scatterers 
become unlikely as the size of each scatterer goes to zero and the Poisson nature of the scatterer 
locations which means that knowing the location of one scatterer does not give information 
about the locations of the other scatterers. Since analyzing the random Lorentz gas directly 
is beyond the capability of current techniques, this random flight model is commonly studied 
in both the mathematics literature (RT991 Vys06| IBFS001IBCLDM021 IBR14| and the physics 
literature |ADBLP10l ldWB041 lvBM051 IMvB04| to gain insight into the behavior of random 
Lorentz gas models. 

We are interested in the regime in which the particle’s velocity is (typically) large. There 
are several natural examples of this setting. The first is an external field which accelerates the 
particle towards infinity. For example, the influence of a constant gravitational field has been 
studied for the random flight process in both a constant density of scatterers [PW791 1RT99) , a 
variable density of scatterers IBR141 , and for the periodic Lorentz gas in two dimensions [CD09I . 
The random flight model in a gravitational field has also been used as a model for a particle 
percolating through a porous medium |WE82| . A second example of a high velocity limit is 
particle motion in a centered, mean zero, isotropic force field |DK091 IKR06) . The soft core 
Poisson Lorentz gas is an example of this situation. A discussion of how much of the past work 
on particles moving in random force fields fits into the high velocity framework can be found in 


IDK09I . 

A unique contribution of this work, is that we study fields that are in general not mean zero, 
unlike prior work on random force fields. For example, the constant field directed towards the 
origin fits into our setting. In such fields, conservation of energy implies that particles have 
bounded trajectories. We are interested in the behavior of the system as the velocity of the 
particle, or equivalently the energy of the system, is sent to infinity. Our primary focus is to 
study the free path of the particle, which is the path of the particle between two reflection times. 
The study of the free path has played a key role in much of the work on the hard core Lorentz gas. 
To highlight a few examples, the free path played an important role in the recent derivation of the 
Boltzmann-Grad limit of the periodic Lorentz gas |MSlla : MSilbl . which is substantially more 
complicated than the corresponding result for the random Lorentz gas. Furthermore, it has also 
played a central role in the derivation of a diffusive limit for the two dimensional periodic Lorentz 
gas in a gravitational held |GD09| . For these reasons, the free path has become an object of 
interest in its own right, see e.g IMS 1 01 IMS14al IMS14bl IWenl21 IBZ071IBGW981IGW001ICG03) . 
We determine precisely the asymptotic behavior of the free path of a particle whose domain 
of motion is predominantly in the high velocity regime. The primary complication that arises 
in our study is that the domain of the particle is determined by the external Held and, as a 
result, the particle’s velocity is not bounded away from zero. If the particle enters a part of the 
domain where it is moving slowly, then there can be a long time before a reflection occurs, and 
we must control this in order to analyze the asymptotic behavior of the mean free path. As 
a consequence of our analysis of the free path, we obtain several diffusion approximations that 
illustrate different aspects of the particle’s behavior. 

We remark that our techniques may be of interest beyond the study of the Boltzmann-Grad 
limit of a Lorentz gas. The random flight process we study falls into a general class of models 
known as transport process. Models similar to the ones we study here are used in many areas 
of physics and chemistry as microscopic models for particle behavior. For example, they arise 
in neuron transport jIAB99I . the study of semiconductors |BAD96] . electron transfer dynamics 
|HW131 IWW05] , and numerous other Helds |BR.0Q1 |Deg98[ IDGL811IDLLS13) lHLB+09] . 
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1.1. The Model. We are primarily interested in the process in dimension three, which is the 
most interesting spatial dimension and also a dimension where certain technical simplifications 
occur. Let us start with a more detailed description of the Boltzmann-Grad limit. Let IA : 
R 3 —> R and g : R 3 — > (0, oo) be differentiable. The function U will serve as the potential for a 
conservative force and g will be the density of scatterers. Fix an energy level E and consider a 
particle moving in the potential IA with total energy E. By conservation of energy 

-+W(y(£)) = -E, 


TO V 


where y is the particle’s position and v is its velocity. In this setting, the Boltzmann-Grad limit 
can be obtained as follows. Assume spherical scatterers with radius 1/R are placed so their 
centers are the points of a Poisson process with intensity R 2 g. Further assume that for almost 
every initial condition the trajectory of a particle moving in the potential U with total energy 
E is open and has infinite length. The arguments of |Spo78[ |Spo88| can easily be adapted to 
include the external field and produce the following result: if the initial position and velocity of 
the particle is absolutely continuous with respect to Lebesgue measure on the constant energy 
surface then the distribution of the position and velocity process of the particle converges as 
R oo, in the sense of convergence of finite dimensional distributions, to a Markovian random 
flight process (X(£), V(£)) t >o with generator 


(1.1) D/(x, v) = v • V x /(x, v) - — S7JA • V v /(x, v) 

m 

+ 2( x )IMI / (/( x , IM|u)-/(x,v))cr(du), 
is 2 

where a is the normalized surface measure on the unit sphere S 2 := {x : ||x|| = 1}, see |Spo88| . 
A similar result holds in dimensions other than 3, but instead of the integral being against the 
normalized surface measure it is against a kernel that depends on v. 


Remark 1.1. The assumption that for almost every initial condition the trajectory of a particle 
with mass to moving in the potential IA with total energy E is open and has infinite length 
is essential to the derivation of the Boltzmann-Grad limit. The intuition for its necessity is 
that it is possible for a periodic orbit to avoid all of the scatterers in a random arrangement. 
If this happens, the effect on the limiting process is that one must keep track of whether or 
not an orbit has closed, which creates long range dependencies in the limiting process. The 
effect of this phenomenon on the Boltzmann-Grad limit has been the subject of study in both 
the mathematics and physics literature |BHHP96l IBMHH971IBHPH0H IDR041IKS98) . In these 
cases one may still obtain a limiting process with generator of the form (ED, but additional 
randomness must be added to array of scatterers. For example, they can be allowed to move 
(but not transfer energy), or they can be arranged according to a space-time Poisson process 
where a scatterer comes into existence, exists for a certain amount of time, and then disappears. 

The process (X(£), V(£))*>o can be constructed in the following way, which explains the name 
“random flight process”. Let (y(x 0 ,v 0 ,£))t>o with x 0 ,v 0 £ R 3 and to||v 0 || 2 /2 + U(x 0 ) = E 
denote the solution to the initial value problem 



( my" = 

-V«(y) 

(1.2) 

{ y(o) = 

x 0 


l y'(o) = 

v 0 - 


We construct our process ((X(£), V(£)), £ > 0) recursively as follows. Set (X(0), V(0)) = (x 0 ,v 0 ) 
and let T 0 = 0. For k > 1, assuming we have defined ((X(£), V(£))) 0 <t<T fc _i, we let U fc _! be 
independent of this part of the path and uniformly distributed on S 2 and let T'k satisfy 


(1.3) P (T k - T fc _! > £ | U fc _ 1 ,((X i ,V t )) 0 < t < Tfc _ 1 ) 

= exp J g (y(X(Tfe_ 1 ), || V'(r fc _ 1 )||U fc _ 1 , a)) ||y'(X(T fc _ 1 ), || V(r fc _!)||U fe _r, s)|| ds^j . 
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For t £ [Tfc_i,Tfc] we then define 

x(t) := y (x(r fc _ 1 ) ! ||v(r fc _ 1 )||u fc _ 1 ,t-r fc _ 1 ) 

(L4) V(t) :=y'(X(r fc _ 1 ) ) ||V(T fc _ 1 )||U fc _ 1 ,t-T fc _ 1 ) 

We remark that under very mild assumptions T*, —> oo a.s. and thus this defines the path of 
the particle for all times. Intuitively, T*, defines the kth. reflection of our particle by a scatterer. 
By studying the free path of the particle, we mean to study, in particular, the conditional law 
of (X(£)) 2 ’ 1 < t < 7 ’ 2 given X(Ti) = x under an appropriate scaling of the parameters in the model. 
Note that the path of the particle on the time interval [Tj, T 2 ] is the path between two consecutive 
reflections. This is why we call it the free path. 

At this point we make a simple but important observation. By conservation of energy, 

l|V(t)|| = ^(i?-«(X(i)), 

so that, if we define 

(1.5) v(x) = J' —{E — W(x)), 

V m 

then 

(1.6) X(t) = y(X(Tfc_ 1 ),«(X(rfc_ 1 ))U fc _ 1 ,t-r fc _ 1 ). 

Since Ufc_i is independent of ((X(f), V(f)))o<t<T fc _u this implies that if we define X& = X(Tfc), 
then (Xfc)fe>i is a Markov chain. That the index in this chain starts at 1 is an artifact of our 
deterministic choice of V(0). If instead of choosing V(0) = Vo in the construction above we took 
V(0) = u(X(0))U, with U uniformly distributed on S 2 , then (X*,) fc>o is a Markov chain and its 
transition operator is 

(1.7) P/(x) = E [/ (y (x, «(x)U, AT(x, U)))], 
where IV (x, u) is a random variable with distribution 

(1.8) P (IV (x, u )>t) = exp J g[ y(x, u(x)u, s)]u[y(x, n(x)u, s)]ds^ , 

and conditional on U = u, 7V(x, U) is distributed like 7V(x, u). 

In order to avoid the notational difficulties and technical assumptions on the geometry of the 
set {x £ R 3 : U(x.) = E} that come from working with a general potential U, we restrict ourselves 
to the spherically symmetric case where W(x) = lA{r) and g(x) = g{r), where r = ||x||. In this 
case angular momentum is conserved, which implies that the motion of the particle between 
collisions with scatterers takes place in a plane. Consequently, we may restrict our attention to R 2 
for some of the analysis of the free path. This lets us work in polar coordinates for these aspects of 
the analysis. Throughout let R + := [0, oo), No := {0,1, 2, ... } = NU{0}, and C k (S) be the set of 
real valued functions on S having k continuous derivatives. We let U{r) be the potential and we 
assume our particle has mass m and total energy E. Let r(t) := ((r(r 0 , 6, t), a(r 0 , on, 0, t)),t > 0) 
denote the trajectory in polar coordinates of a particle with mass m > 0, total energy E, initial 
position (ro,ao) and whose velocity vector at time 0 makes angle 9 with the radial vector. 
Consequently r(f) is the solution to the equations of motion in polar coordinates for a particle 
of mass m in the potential IA : 

. . 1 d , 9 .. ... . 9 . d r U(r) 

(1.9) - — (r 2 a)e a + (r — ra 2 )e r =- -e r 

r at m 

with initial conditions 

r (0) = (ro,ao), 
r(0) = r’(ro) cos 9, 

d(0) = ^2 sind, 
ro 


( 1 . 10 ) 
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where e r is the radial unit vector, e a is the angular unit vector and 6 is the angle the initial 
velocity makes with the line segment from the origin to the point r(0). We remark that in this 
setup it is generally possible for the trajectory to hit the origin in finite time, at which point 
the equations of motion are no longer well defined. In this case we continue the path in such a 
way that momentum is continuous, if possible and such that angular momentum is conserved. 
The extension is most easily described in Euclidean coordinates, since it may involve the angular 
component reentering from ±oo in polar coordinates. Suppose that y(t) is a trajectory that hits 
the origin for the first time at time T > 0. We extend the path to the interval [T, 2 T] by defining 
y (t) = —y(2 T — t) for t £ [T,2T\. Since ||y(2T)|| = ||y(0)|| ^ 0, we can continue the path by 
the equations of motion until it hits the origin again, at which point we repeat the extension 
procedure. A straightforward calculation using the radial symmetry of the potential shows that 

^(-y(2 T - t)) = T - t)), 

so our extension solves CE3 away from the origin and it is easy to see that y' is continuous at 
T if lim t^Ty'(t) exists. In this way, we guarantee that all initial conditions lead to trajectories 
that exist for all time. 

Let v(t) := v(r(t)) denote the speed of the particle as a function of time. By conservation of 
energy, we have 

( 1 . 11 ) ^M +ll (r) = E 

and, as a result, 

(1.12) v(r) = \J -|( E-U{r )). 

We can also write v(t) as a function of angular and radial velocity: 

(1.13) v(t) = \J r 2 {t) + r 2 (t)a 2 (t). 

1.2. Heuristics. In this section, we give heuristic arguments to motivate our scaling. Consider 
the following simple situation that encapsulates the phenomena in which we are interested. In 
particular, suppose that m = 2, U(x) = ||x||, so that IA produces a uniform acceleration towards 
the origin, and that g(x) = 1. Suppose we start the particle at the origin with speed vo = 1 
moving in direction g. By conservation of energy the total energy of the particle will be 

E(vq) = mvl/2 = 1. 

Since IA produces a uniform acceleration towards the origin, the maximum distance the particle 
can move away from the origin is 1. Consequently, there is no way to rescale the trajectory of the 
particle to obtain a diffusive limit: Any scaling of space will cause the trajectory to degenerate, 
and scaling time alone will produce jumps. However, we can obtain a diffusive limit if we let 
vo (and, consequently, E ) go to infinity. This is the type of high velocity limit considered in 
[DK091 KR06 :. In those papers, however, the high velocity was used to prevent the particle 
from being trapped in a bounded domain, while in our setting the particle is still trapped in a 
bounded domain, but the size of the domain goes to infinity as vq does. 

Treating vo as a parameter, the energy becomes a function of vo, specifically E(v o) = Vq is the 
maximum distance the particle can travel from the origin. Let (X^° (i))t>o be the trajectory of 
the random flight process with initial conditions X'" 0 (0) = 0 and initial velocity X 1 ’ 0 (0) = i>oU, 
where U is uniformly distributed on S 2 . In order to obtain a diffusive limit for (X' u °(t)) t > q, we 
rescale space so that its maximum distance from the origin remains constant. In particular, we 
look at (vq 2 X Vo ( t)) t > o, which travels at most distance 1 from the origin. This is the process for 
which we analyze the behavior of the free paths. 

Since the velocity of the particle is going to infinity, the time between reflections is going to 
0, so we can expect a diffusive limit if we scale time appropriately. Since g = 1, the distance 
(X v °(t))t>o travels between reflections has an exponential distribution with rate 1, so the dis¬ 
tance (i> ( j' 2 X ,;o (t)) t >o travels between reflections is of order Vq 2 . Writing down the equations 
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of motion for (vq 2 X. V0 (t))t >o between reflections, we find that its trajectories are the same as 
those governed by the potential U v °{x) = Vg 2 U(x) moving with total energy E v ° = Vq 2 , initial 
velocity Vq , and the scattering rate as measured along these paths is g v °{x) = Vg. 

By Brownian scaling, in order to obtain a diffusive limit, we expect that we should scale time 
so that there are approximately Vg reflections per unit time. Typically, we expect the particle to 
be distance of order Vg from the origin, and if X”° (t) = VgX and t is a reflection time, the time 
until the next reflection is approximately 


1 


1 _ 1 
V v o -^oll x ll W 1 - ll x ll 


Thus the amount of time for Vg reflections to occur is of the order Vg. This suggests looking for 
a diffusive limit of (v^ _2 X L ’° (vQi))t>o as vg tends to infinity, which is what we undertake. 


1.3. General scaling and assumptions. In general, we allow g to be a small perturbation of 
the constant density and IA to produce a small perturbation of the constant field, small being 
relative to the speed of the particle. In particular, we will assume the following conditions. 

(Al) The path of the particle (X"(f)) t >o evolves as the position of a random flight process 
where the density of scatters is 

g n ■■= Vng , 

the potential energy of the field 


and the total energy of the particle 


U n := -±=U 

\/n 

E n := —r=E, 
yjn 


where g, U , and E are fixed and independent of n and g(x) = g(||x||) and U(x) = W(||x||). That 
is, the density of scatterers and potential are spherically symmetric. By (11.121) we can write the 
speed as a function of E and U ; consequently the speed v is also rescaled as 

1 


:= 


d/4 


){r). 


The trajectory of the particle with these rescaled parameters we denote by 

(yn(x,u„(x)u,t)) t >o 

when we are working in with the path in 3 three dimensions and by 

((r n {r o ,O,t),a n (ro,a o ,0,t)),t > 0). 

when we are considering the path in polar coordinates with respect to its plane of motion. When 
no confusion will arise, we will leave the dependence on initial conditions implicit. 

In general, subscripts or superscripts of n refer to distributions relative to these scaled param¬ 
eters. For example, corresponding to m is the random variable N^ (x, u) with distribution 

r t 


(1.14) P (A (n) (x, u) > tj =exp J ff n [y n (x, u n (x)u, s)]w n [y n (x, u n (x)u,; 
and corresponding to (11.71) is the Markov chain (X^)fc> 0 with transition operator 


0 )di 


(1.15) 


P„/(x) = E [/ (y„ (x,»„(x)U,ffW(x,U 


where U is uniformly distributed on S 2 and conditionally given U = u, N( n ' > (x, U) is distributed 
like N^ n \x, u). Similarly, (X n (f)) t >o will be constructed as in (11.61) we the rescaled parameters. 


Remark 1.2. To make the connection with are heuristic arguments in Section 11.21 if we take 
vo = n 1 / 4 , g = 1 and U(x) = ||x|| then (X”(f)) t > 0 (vq 2 X. Vo (t)) t >o- Thus the scaling we do 
here accounts for increasing the initial speed and scaling space to keep the particle’s trajectory 
contained. This is the scaling under which we analyze the free path of the particle. We will scale 
time separately when we consider diffusive limits. 
















FREE PATH IN AN EXTERNAL FIELD 


7 


(A2) We assume that for all n £ N the process (||X"(t)||, t > 0) evolves in a domain V C R+ 
where T> = K+ or T> = [h-,h+\, where 0 < h- < h+ < oo. The domain T> is chosen so that 
E — U(r) >0 for all r £ T>°. This is equivalent to v(r) > 0 for all r £ T>°. 

We require the particle to have positive speed in the interior T>° so the time between reflections 
approaches 0 as n goes to infinity. This way we obtain a nontrivial diffusion limit. 

(A3) On the boundary dV we have the following assumptions. 

• If fc_ > 0 then 

U(h- 
— d r lA(h- 

• If h + < oo then 

U{h+) = E 
-drU{h+) < 0. 

The conditions on IA force the speed of the particle to be zero at the respective endpoints /i_ 
and h+, while the conditions on ^ ensure the force field at the endpoints is pointing towards 
the interior T>°. Having zero speed at the boundaries and having the force field point ‘inwards ’ 
prevents the particle from leaving the domain T>. 

• If h- = 0 then 

E-U{ 0) > 0. 

This condition ensures the particle is not trapped at the origin. 

• If h + = oo we require that for any e > 0, infr h_+e,oo){E —U[r )) > 0. 

This condition implies that for any e > 0, inf[/ l __|_ €jtx) ) v(r) > 0 which shows that the reflection 
rate does not go to 0 as the process goes to infinity. 

(A4 )U&C 1 (V) 

This smoothness assumption ensures that the velocity and the acceleration of the particle 
depend continuously on the position r(t). If T> = [0,h] we can relax this condition to IA £ 
C 1 ((0, h\) so that we can allow potentials of the form lA(r) = — ^ which are not defined at 0. 
(A5) The density g is spherically symmetric and satisfies 

g eC(v)nc\v°) 

together with 

inf g(r) > 0. 
rev y 

We require some smoothness from g because the diffusion limit we get depends on the derivative 
g'. The second assumption is needed because we do not want to have regions where the the 
reflection rate goes to 0. If one allows g to approach 0, then the free path may scale differently 
when the process is started in these reqions. Examples of such situations are dealt with in Section 
8 of lBRl4l . 


>-V»o 


1.4. Results. Our first result considers the asymptotic behavior of the steps of the free path 
Markov chain ((X£, Tjf))k>o, where X£ is the location of the particle at the time of the fc’th 
reflection and, anticipating our diffusion approximations, is the time of the fc’th collision. 

This Markov chain has transition operator 


(1.16) 


Qn/( x , t) = E 


/ (y [x, «(x)U, iV"(x, U)], t + n- 3 / 4 A"(x, U)) 


The following result characterizes the asymptotic mean and covariance structure of the free path 
chain. 


Theorem 1.3. Let 


Mx,t) := nE[(X?,7?) - (x,i) | (X£,T 0 n ) = (x,i)] 
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be the scaled drift of (X£, Tff)k>o, and let S be a compact subset ofV°. Then 


lim sup 

71-5,00 (x,i):||x||eS,i>0 




-1 


2W(x) , V.g(x) 

ff( x ) 


s(xMx) 


^3g(x ) 2 mu(x ) 2 
uniformly on {x : ||x|| £ S'} x [0, oo). Furthermore, let 

<y(M) “ nE [(Xft - *0 (X^- - | (Xg, Tf ) = (x,i)] 


= 0 , 


<«(x,t) := nE [(Xft - a*) (71" - i) | (XJ,7?) = (x,t)] 


Then 


o£ it (x,t) := nE 



i) 2 |(X^,V) 


(x,i) 


( 1 ) 

( 2 ) 

(3) 


2 

<W x .O- 3^)2 

lim sup o’L t (x,t)= 0 . 

71-700 (x,t):||x||eS,t>0 

lim sup u„ f (x.,t) = 0 . 

71-700 (x,t):||x||eS,t>0 ’ 


lim sup 

71-700 (x,t):||x||es,t>0 


= 0 , where Sij 


t{i=j}. 


This theorem is a combination of Lemmas 12.101 and 12.111 below. Heuristically, this says that 
both U n and g n impart a drift towards areas where the corresponding function has a smaller 
value. That is, the particle prefers to move towards areas where it travels quickly and where there 
are few scatterers, which may be competing influences. It is interesting to note that the potential 
U does not appear in the covariance terms. It can, and will, effect the diffusion coefficient of the 
limiting diffusion only through an overall time change. 

With Theorem EH in hand, it is straightforward to obtain the following diffusion approxi¬ 
mation, which separates out the effects of the particle’s position at times of reflection from the 
effects of speed at which the particle is moving. Let = {x £ R 3 : ||x|| £ T>}. 


Theorem 1.4. Let (X t )t> o be a diffusion onV 3 whose generator G acts on functions f £ C 2 (fDf) 
with compact support in T >3 by 


(1.17) 


G/(x) 


1 

3ff 2 ( x ) 


A/(x) 


1 ( Vff(x) 

3<? 2 ( x ) V s(x) 


2 V77(x) 
rov(x ) 2 


V/(x) 


and killed if/when 77 := ||X|| hits the boundary ofT>. Consider any l,u £ T>° with l < u and 
start the process ((X£,7)."),fc £ No) at (x, 0), where l < ||x|| < u. Define the stopping times 


r ”« := inf{fc£ No :|ra£ M} 


and 

n,u := infji > 0 : \\X t \\ [l,u]}. 

Then, as n —> 00 , the family of continuous time processes ((X” nt j At „ , 7^ t j Ar n ),t > 0) con¬ 
verges in distribution on the Skorokhod space 7?(M+,K) to the diffusion 


with initial position (x, 0 ). 
times, n~ 1 rjf u =>• n iU . 




ftATi' U 


ds 


tAT iiU j 


,t> 0 


0 g(X.)v{X a ), 

This convergence happens jointly with the convergence of the hitting 


This theorem is proved in Section [3l where it is restated for convenience. 


Remark 1.5. The convergence above also holds without stopping near a boundary point of D 
that is inaccessible for the diffusion 77 = ||X||. See Section [5] for how one can determine when a 
point is inaccessible. In particular, the boundaries are inaccessible for the constant acceleration 
towards the origin and for Newtonian gravity centered at the origin. 
























FREE PATH IN AN EXTERNAL FIELD 


9 


Remark 1.6. This type of limit theorem, that gives joint convergence of a process observed when 
it changes direction and the time spent between such changes up to that point, is typical of the 
continuous time random walk literature, see e.g. | MS14cl . Although our models can naturally 
be viewed as transport processes, our results on the free path of the particle suggest an approach 
to diffusive limits that has more in common with that used to study Continuous Time Random 
Walks (CTRW) [MS041 IMS14c| than that typically used to study transport processes. CTRW 
models also arise in many physical applications (see e.g. [WZL + 14l IMK001IBCDS061ISBMB03 ] ) 
and have been particularly useful in studying anomalous diffusions. The advantage of this 
approach is that, on the diffusive scale, we may easily distinguish between the effects of the 
particle’s displacement between collisions and the effects of the speed at which the particle is 
traveling. In this way, our approach is similar to the one taken in jBR14| . and our work can be 
seen as an extension of the methods used there. This approach is also similar to the approach 
in [CD091 to studying the two dimensional periodic Lorentz gas in a gravitational field, though, 
of course, our setting is far simpler. 

By inverting the time process {Tff)k>o and showing that the free path is typically not far from 
the straight line between its endpoints, we can arrive at a diffusion approximation for (X"(f))j>o, 
the trajectory of the particle. In particular, we obtain the following result. 


Theorem 1.7. Let (X”(t)) denote the trajectory of the particle and fix l,u £ T>° with l < u. 
Suppose that ||X”(0)|| £ ( l,u ) and define 


il u := inf{t:||X"(i)||£M}. 

We have the following convergence in distribution on the Skorokhod space U(R+,] 

(x”((n 3 / 4 f) A {X{m A n,u)) t 

where X is as in Theorem no and fl is the time change given by 

ds 


n> o 


fl(t) :=1 


(*) 


/0 g(X(s A Tl,u))v(X(s A Tl^u)) / 
and I is the inverse operator defined by I{f)(t) = inf{s : /(s) > t}. Furthermore, the time 
changed process (T(f l(t)),t > 0 )) is a diffusion process whose generator Sf acts on functions 
f £ C 2 (T> 3 ) with compact support in V 3 by 


»/(*)=«(*)»WG/W = ^A/(x)- 


3 s(x) 


35 (x) V ff(x) mu(x) 2 


V/(x) 


This theorem is proved in Section [H where it is restated for convenience. 


Remark 1.8. The stopping at u and/or at l can be removed when the left and/or right boundary 
points of V are inaccessible for the diffusion 1Z. See Remark [L5l and Section [5] 


Remark 1.9. The cutoffs are necessary because of the generality the potentials and scattering 
densities we allow permit very different behaviors at the boundary. In some cases we the bound¬ 
aries to be inaccessible (see Section [5]), while in other cases we expect the boundaries to be 
reflecting as in jRT991 and in yet others, like U(r) = — r -2 , the origin should trap the particle. 
We leave a detailed investigation of the boundary as an open problem. 

Remark 1.10. Our results have much in common with |Cos91| . but there are key differences 
at both heuristic and technical levels and, as a result, there is no overlap between the precise 
results. For instance, our models do not satisfy the underlying technical assumptions of ICos91j . 
The function Q -1 from |Cos91| . which controls the reflection rate, has to be differentiable up to 
the boundary of the region where the motion takes place (see assumption (H2) from |Cos91j ). 
In our case we have singularities at the origin and when T> = [0 ,h\ we also have singularities 
at the boundary of the domain. Furthermore, the speed of the particle must be bounded away 
from 0 and 00 in |Cos91| . while there is no such restriction in our model. This is related to the 
fact that [Cos91| restricts the particle’s motion by imposing a reflecting boundary while in our 
case the domain of the particle’s motion is completely determined by IA and E. 
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Remark 1.11. The form of the generator EU puts the random flight process in the class of 
processes known as transport processes. There is a substantial literature on diffusion approxi¬ 
mations of transport processes and, more generally random evolutions. Very general approaches 
and results can be found, for example in [ Kur73l |Pap75| [ BLP79 i Cos91] . In order to apply these 
general techniques the generator needs to be smoothed to control for regions where the parti¬ 
cle’s velocity is low. Nonetheless, it is possible to adapt these approaches to obtain Theorem 
EH However, our results on the free path are stronger than the diffusion approximation results 
and cannot be recovered by these general methods. Moreover, these methods do not clarify the 
connection with the free path chain. 


2. The free path of the particle 


In this section we analyze the jumps of the free flight chain ((XjJ, T£))k>o whose transition 
operator is given by (11.161) . Our proofs are complicated by the fact that for a general potential 
U it is not possible to explicitly find the trajectory of the particle in the absence of collisions and 
this makes it difficult to determine how long this path spends in domains where it is traveling 
slowly. Our first lemma will give us local control over the how the radial part of the path behaves. 


Lemma 2.1. For any ro £ TV and 6 £ [ — 7r, 7t] 

r(t) := r{r 0 ,9 , t) = r 0 + v 0 cos(6») • t + f dr^(K r )) 

2 \ m 


vl r l sin 2 (9) \ 2 
r(r ) 3 ) 


for some 0 < t < t depending on ro, 0 and t. 

Proof. By definition r (t) := (r(ro, 9, t), a(ro, ao, 9, t)) is the solution to the equations of motion 
in polar coordinates for a particle of mass rn in the potential U. Equation (11.91) implies 

(2.1) --§(r 2 d)=0 

r at 

and 


( 2 . 2 ) 

As a result of (12.11) and (11.101) 

(2.3) 

By (12.21) and (12.31) 

(2.4) 


f — rd 2 = — 


d r U (r) 


d(t) = 


v(ro)ro sin0 


r{t) 2 


Ht) = 


-~d r U(r(t)) v(r 0 ) 2 rQ sim(0) 


m r(f) 3 

Taylor expanding r(t) and using (12.21) . (11.101) together with (12.31) yields 


(2.5) 


, , 1 /— d r U(r(r)) ri(ro) 2 r-o sin 2 (0) 

r(t) =r 0 + v{r 0 ) cos(0) • t + - V K U + V J ° 3 V ; 

2 \ m r{T)° 


t 2 . 


□ 


Remark 2.2. Throughout the remainder of this section we let S, S' be closed intervals satisfying 

S C (S")° C S’ C V°. 

Also, for any <5 >0, we define 

(2.6) A s(S)= inf inf {t: ||y„( x , u, t) - x|| > 6} 

(x,u):||x||eS,uGS 2 t> 0 

to be the shortest time it takes for the particle’s displacement from its initial position to be at 
least 8 when started inside S. 

If 

0 < 8 < d(S, ( S ') c ) := inf {\x — y\ : x £ S,y £V\ S' 7 } , 
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then for all (x, u, t) £ {x : ||x|| £ S} x S 2 x [0,A^(,S')] one has 

||y„(x,u,f)|| £ S'. 


Since ||y„|| = n 1//4 u(y), we have 


Ag(S) > n 1/4 ■ 


su Pr 0 eS’ v ( r o) 


> 0 


where sup rog s/ u(r’o) < oo since S' is bounded away from dD. 

The next lemma shows that the scaled trajectories converge to the starting point, uniformly 
over a fixed time interval. 

Lemma 2.3. Fix T > 0. Then 

lim sup ||y„(x, u, t) - x|| = 0 

n ^°° (x.u,t)6{x:||x||eS}xS 2 x[0,T] 


lim sup 

n_>0 ° (x,u,t)G{x:||x|| GS} xS 2 x [0,T] 


V^y«(x, u, t) 


VW(x) 

m 


= 0 


and 

(2.7) 


lim sup | r n (r 0 , 9, t) - r 0 1 = 0 

n_>0 ° (r-o,e,t)eSx[-7r,ir]x[0,T] 


lim sup 

n_>0 ° (r o ,9,t)eSx[-7r,ir]x[0,T] 


Vnr n (r o ,0,t) - 


d r U(r 0 ) 

TO 


v 2 (ro) sin 2 (0) 


r o 


= 0. 


Proof. The claims about r n are slightly more complicated, so we prove those and leave the more 
obvious claims about y n to the reader. Let <5 > 0 be such that <5 < d(S , ( S ') c ). By Remark 
[U there exists M £ N large enough such that A %(S) > T whenever n > M. Equivalently, for 
n> M we have 

\r n {ro,0,t) - r 0 | < <5 

for all (ro,6,t) £ S x [—7r, to] x [0,T]. This proves the uniform convergence of r n (r 0l 9,t) to r 0 . 
Define 


( 2 . 8 ) 


ip(r o ,r o ,0) := 


d r U(r 0 ) | v 2 (r 0 )r 2 sin 2 (0) 

I Q 

m rg 


and note by (12.411 that 

r„(r o ,0,t) = —if(r n (t),r o ,0). 

\Jn 

S’ is bounded away from 0 and U £ C 1 (D°) imply that if is uniformly continuous on S' x S x 
[—7r,7r] . By construction r n (t) £ S' for all n > M and t £ [0,T]. Because r n (r 0 ,9,t) converges 
uniformly to tq on compact sets, we have 


(2.9) lim sup {\Vnr n (r 0 ,9,t) - V>(r 0 ,r 0 ,6»)|} = 0 

n ^°° (r o ,e,t)eSx[-7r,7r]x[0,T] 


where 

ip(r Q ,r Q ,0) = 

by (12.81) . This completes the proof. 


d r U(r 0 ) v 2 (r 0 ) sin 2 (9) 

TO To 


□ 


Since the free flight chain tracks the process at reflection times, we want to apply the estimates 
for (12.51) between reflections. By (11.31) and the rescaling, we know that for every k, if at Tk the 
particle is at position (r 0 , cx 0 ) and reflects in the direction given by the angle 9, then TJf +l — Tjf 
is distributed like the random variable (r o, 9) which we define by 

P(N ( - n \r 0 ,9) >t) = exp J n 1/4 X(r n (r 0 ,9,s)) ds 


(2.10) 
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where 

(2.11) A (t 0 ) := g{r 0 )v(r 0 ). 

Throughout, we will often suppress the r 0 and 9 dependencies of N' n \ro,6) and write N^ 
or N^ n \9) when no confusion will arise. To go back and forth between Euclidean and polar 
coordinates, note that we have the identity in distribution 

7V (n) (x,u) = d 7V (n) ^||x|| , cos -1 » 

where u) is defined as in (11.81) with the appropriately rescaled parameters. Consequently, 

for studying the asymptotic duration of the free path, there is no loss in studying N^ (r, 6) in 
a uniform way over r and 9. 

For many of our proofs we require estimates that show the time between reflections approaches 
0 with high probability as the scaling factor n goes to infinity. This will allow us to apply the 
local estimates from the expansion of r(t) in Lemma HOI and to show the free flight process does 
not undergo large jumps. We first prove some bounds on the moments of N^ n \r 0 ,9). 


Lemma 2.4. The family 

{n 1 / 4 -/V’(”)(ro, 9) : (ro,9, n) £ S x [— n, n] x N} 
is bounded in L p for 1 < p < oo. 


Proof. We will assume that T> = [0, h] or T> = [0, oo). The cases V = [h-, h+] and [h, oo) can be 
treated similarly. 

Case I: T> = [0, h] 

By Assumptions (A3) and (A4) there exist S > 0 and Us > 0 such that 

(2.12) min \d r U(r 0 )\ > Us- 

[h—8,h] 

Let S C D° be a compact set, let rj > 0 and assume our particle enters the annulus A v with 
inner radius h — g and outer radius h at time to ■ By (11.121) 


v (h — Tj) = \ —(E-U(h-rj)). 
m 


f(t> - JJPAll+ r(t) # w < JJCAll + 


(2.13) 

Using (12.21) and (12.31) 

(2.14) 

Let 

t e := inf{s > to : r(s) = h — 77 } 

be the time when the particle exits the annulus. By Assumption (A3) 

U(h) = E. 


m 


TO 


Since U is continuous this means that we can make v(r) as small as we like if we are close enough 
to r = h. This together with 02.121) . (|2.13L and 02.141) implies that there exist 7 > 0,m 7 > 0 
such that 

r(t) < —to 7 

whenever r(t) £ [h — 7 , h\. Set 77 = 7 . Note that t e — to is bounded above by the time it would 
take a particle started at r = h — 7 with speed r(fo) > 0 pointed along the radius and with 
acceleration r = —to 7 < 0 to return to r = h — 7 . Thus, 

(2.15) t e - t 0 < 2—- L —^ < 00 . 

to 7 


Next, define 

t r := inf{s > t e : r(s) = h — 7 }. 

This is the first return time to the annulus A y . We want to bound t r — t e below. If r(t 0 ) < 0 
then the particle would not spend any time in the annulus Ay, that is t e — t 0 = 0. Therefore, we 
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can assume that r(to) > 0. By conservation of angular momentum and conservation of energy 
we have that v(to) = v(t e ) and a{t e ) = a{to). This implies that f(t e ) = —r(to) < 0. Since f is 
finite on S we immediately get that 


(2.16) 


t'p tg ^ 


2 r(t 0 ) 


> 


2 r(t 0 ) 


> 


2 r(t 0 ) 


sup s |f| - gHPsMlll + supg |rd 2 | “ sup s \0M(r)\ + supg v^r) 
Combining (12.151) and (12.161) . 

o U*o) 


> 0 . 


/r, 1 rj\ t e ^0 

(2.17) sup -= sup 

Sx [— 7T,7r] tr t e Sx[— 7T,7r] 


2r(tp) 


aupg \d r U(r)\ | =U p^ " 2 Q) 


sup s \d r U(r)\ , „„„ V (r) 

< - +SU P^^ <0o 

m 7 


where we assume that if to = oo then t e — to = 0 and 

te t 0 _ ^ 

t r t e 

Note that the upper bound from (12.171) is the same if we use the rescaled trajectory of the 
particle (r n (t), a n (t)). This follows from the way we found m\ and by Assumption (Al). 

By Assumption (A2) we know that 

inf v(r o) > 0. 

Suppose that t > sup s (t e — to) + infs(t r —t e ). 

Using (12.171) together with the fact that the worst case scenario is when the particle spends the 
longest possible time in the ‘bad region’ A 7 and the least amount of time in the ‘good region’ 
T>\A 1 , we have for some c > 0 that 

,-t 


(2.18) 


P{ 7 V(™) > t} = exp J n}^g(r n )v{r n ) ds 
( ri 


< exp — 


n 1 / 4 inf g(r) inf v{r) ds 


= exp - 


-tn 1 / 4 in f q(r) inf v(r) 
1 + c S w 3 


which decays exponentially in n as n —> oo as long as t is large. Therefore, 

p . •“ 


E 


(n 1/4 M n) y 


(2.19) 


= J pt^ 1 P (n 1/4 M n) > i) dt 
< J pt p ~ x exp tinf g{jr) inf v(r)^j dt 


< oo. 


Since the bound above does not depend on r, 9 or n we get that 


sup E 

S X [ —7T,7r] XN 


^ 7l 1 / 4 ^\r( ra )( r °. e ))y 


< oo. 


Case II: V = [0, oo) 

Set g m i n := inf re x> g( r )- We know by Assumption (A3) that there exists 5 > 0 such that 
inf-p v(r) > 6 so by (12.101) 

P{AT (n) > t} < exp (-n 1/4 6tg min ^ 

which, like before, forces 


E 


(n 1/4 lV(")) ? 


= J ptP- 1 P (n 1/4 N^ > tj dt 
/»00 

< / pt p ~ l exp (—6tg m i n ) dt 
Jo 


Jo 

< oo. 


( 2 . 20 ) 
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Similarly to Case I, because the bound above does not depend on r, 9 or n we get that 

r / , w rm\Pl 

sup E 

SX [ — 7T,7r] XN 


(n 1/4 iV^ r °’ e »y 


< oo. 


□ 


Lemma 12.41 also provides us with the following corollary which will prove useful in showing 
that the probability that (ro, 9) is larger than any fixed value decays rapidly as the scaling 
parameter n goes to infinity. 

Corollary 2.5. For all k £ R and for all e > 0, 

lim sup n fc P \N^ n \r 0 , 9) > e| = 0 

n_> ' 00 (r 0 ,6>)eSx[-7i-,7r] > 

Proof. This is immediate from Lemma 12.41 

□ 

In addition, we have also the following corollary that shows that tails of the moments of 
N^ n \ro,9) decay rapidly as well. 

Corollary 2.6. For all k € R, Z > 1 and e > 0 


lim sup n k E 

n—loo (r 0 ,e)eSx[—7r,7r] 


(N n (r 0 ,9)) l l {NW>E} 


= 0 . 


Proof. Cauchy-Schwarz together with Lemma LTTl and Corollary 12.51 yield the desired result. 

From Lemma m we have the following estimates on the moments of . 

Lemma 2.7. Let 1 < p < oo. Then 


□ 


(n p/ 4 [lV (n) (ro,#)]^ - f pt p 1 exp (—g(r 0 )v(ro)t) 


dt 


= 0 . 


(2.21) lim sup E 

ra-loo ( rOi 0)eS x [-7r,7r] 

Proof. First let M £ R + and note that the truncated moments E [n p ^ 4 N^ n \ro,9) p ] A M p ] can 
be written as 


E 


( 2 . 22 ) 


n J>/ 4 /v(")( rO) 0 )P a M p ] = p t p _1 P \n 1/4 N^ n \r 0 ,9) AM >t 


dt 


r M 


= p / tP- 1 P n 1/4 iV (n) (r o ,60 > t 


dt. 


Making the change of variables u = n 1 // 4 s, we have 

(2.23) P ^n 1 ' / 4 iV( ri )(ro, 9) > t^j = exp J g(r(r 0 ,9, u/n 1 // 4 ))n(ro, 9, u/n 1 ^ 4 )^ du. 

Both r = r(ro,9, t) v = v{tq, 9, t) are continuous functions on O := S x [— n, 7 r] x [0, M}. Since O 
is compact, r and v are in fact uniformly continuous on S. By Assumption (A5) g is continuous 
on S, and therefore is uniformly continuous. This implies that g o r is uniformly continuous on 
O. By Lemma [?~3l it follows that 

(2.24) lim g(r n (r 0 ,6,u/n 1/4 )) = s(r(r o ,0,O)) = g(r 0 ) 


and 

(2.25) 


lim v(ro, 9, u/n 1 ^ 4 ) — v(ro, 9, 0) := v(ro) 


both uniformly on O. As a result, P ( n 4 ^ 4 N^ (ro, 9) > f) converges uniformly to exp (—g(ro)v(ro)t) 
on O, which implies 


(2.26) lim sup 







E 

n p/4 N^ n \r 0 ,9) p A M p 

- / pt p_1 exp (—s(r 0 )n(r 0 )f) 

7T,7r] 



Jo 


= 0 . 
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It is then standard to extend the result to the expectation without the truncation. 


□ 


For the cases p = 1,2 we have by Lemma 12771 


(2.27) 

and 

(2.28) 


lim sup 

n —,►oo / o \ c r 


lim sup 

(r 0i e)6Sx[-7r,7r] 



E 

n i/ 4 /V(”)(r 0 ,6l) 

1 

7T,7r] 



g{ro)v(r 0 ) 


= 0 


E 


i 1/2 (n< n \r 0 ,9)y 


5 2 (r 0 )u 2 (r 0 ) 


= 0. 


The next result shows that on the event {N^ > e} the fcth moment of the displacement of 
the particle at time 0 and at the first reflection N decays faster than —^ as n —> oo. 


Lemma 2.8. Fix e > 0 , m, k £ N. Then 


lim sup n m E 

" _> °° {(x 0 ,u):||x 0 ||eS,ueS 2 } 


(V n >(u)) 


-x 0 


l{ATW>e} 


= 0. 


Proof. Since ||y'|| = v n (y) we have 


E 


fc if \ K 

(lY (n) (u))-x 0 l {J v(")> e } < E [J v n {y(s))ds\ l {JV W>e} 


< n- k/2 g~t E 


rN^ 


2k' 


9 n(y{s))v n (y{s))ds 


1/2 


P (N^ > e) 1/2 . 


The lemma now follows from Corollary 12.51 and the fact that 

/ 9 n(y(s))v n (y(s))ds 
Jo 

is distributed like an Exp(l) random variable and, consequently, has finite moments. 
Lemma 2.9. Let U be uniformly distributed on S 2 . Then 


□ 


(2.29) 


lim sup 
7l ^ 00 x:||x||e5 


i 3 / 4 E 


7V (n) (x,U)U 


1 / 7»(x) Vg( x ) 

3g(x) 2 v(x) \mr(x) 2 g(x) 


= 0. 


Proof. First, for notational convenience, define the auxiliary function 


(2.30) 


F n (x,u,i):= / g n (y n (x, u, s))v n (y n (x, u, s))ds 
Jo 

= / 5n( r n(||x||, 0, s))u n (r„([|x||, 0, s))ds, 


where 6 = cos _1 (u • x/||x||). In the usual way, we will often suppress the x and u dependencies 
of F n and write F n (t) when no confusion will arise. As previously noted, F„(x, u, 7V( n )(x, u)) 
has an Exp(l) distribution, so 

(2.31) l = E[F n (x,u,ATW(x,u))] 
and 

(2.32) 2 = E[F n (x,u,ATW(x,u)) 2 ]. 

In order to prove uniform convergence, we need to work on a compact set, so we fix a time 
T > 0 and split the expectation and use the fact that E(U) = (0,0, 0) 

(2.33) 0 = E[t/ 1 F„(x,U,7VW(U))l {JvW(u)<T} ]+E[C/ 1 F„(x,U,fV(™)(u))l {JV( „ )(u) > T} ] 
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By symmetry, it is enough to consider the first coordinate. Taylor expanding F n (x 1 u,t) about 
t = 0 yields 


1 ^ 


F n (x, u, t) = g n {x)v n (r 0 )t + ^F„(x, u, r(t))t 2 


(2.34) 

for some r(t) £ [0, t\. Here F n denotes the second derivative with respect to time t. Consequently 


E 


(2.35) 


UiF n (x, U, iV(")(U))l {JV (n) (u)< t } 

= g n (x)u„(x)E [e/ 1 JVW(U)1 {jv( », (u)<t} 
If we substitute this into (12.3311 we can write 


+ -E 




E 


I/ 1 fV(")(U)l {iv( n, (U)<T} 


-1 


(2.36) 


and so 
(2.37) 

n 3/4 E 


g n (x)v n (x) 
1 


E 


UiF n (x, U, fV^(U))l { jvC«) (u) >T } 


-E 


U)) 2 1 {jv( „ )(u)<t} 


C7i7V( n )(U) 


-n 1 / 2 

g(x)v(x 

+ l E 


E 


U\F n (x, U, (U))l {J y(n) (U)>T} 


H 1 F n (r)(iVW( U)) 2 l {ivW(u)<T} 


+ n 3 / 4 E \UiNW{TJ)1 {nM(U) 2T} 


Next we compute the limit as n —> oo of each term in the expansion (12.3711 . The first and last 
term on the right hand side can easily be seen to go to 0 as follows: By applying Corollarv l2.61 
it follows that 


(2.38) 


lim sup n 3/4 E 7V (n) (x, u)l j N (. n ), u)>T} 

n ^°°(x.u):||x||eS L 


= 0. 


Similarly, by Cauchy-Schwarz, (12.3211 and by an application of Corollary 12.51 with k = 1 we have 


n 1/2 E [F n (x,u,JVW(u))l {JV( », (u) > r} J < n 1 / 2 ^E^x,u, NW(u)]F{NW(u) > T} 

( 2 ' 39 ) = v^Y / nP{iV (ra) (u) > T} 

->• 0 

uniformly for (x, u) such that ||x|| £ S as n —> oo. 

We now deal with the middle term. Differentiating equation (12.301) twice yields 

F„(t) = v n (y n (t))Vg n (y n (t)) ■ y n (t) + s„(yn(i))Vu n (y n (t)) • y n (t) 

= n 1/4 y n (f) • (v n (y n (t))S7g n (y n (t)) + g n (y n (t))Vv n (y n 

To evaluate lim rl _ ) . 00 n 1 / 4 r„(t), note that by (11.2D and Lemma [TTiil 


\\n 1/A y n (t) - v(x)u|| = n 1/4 ||y„(f) - y n 

<n 1/A [ ||y„(s)||ds 
J 0 

1 

174 / Vn\\yn{s)\\ds 
1 ' Jo 


n x 
■ 0 


(2.40) 
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uniformly for (x, u, t) £ {x : ||x|| £ S} x S 2 x [0,T] as n —> oo. Differentiating (11.51) shows 

Vt x) =- 

mv(x) 

which toeether with (12.401) . Lemma 12.31 and the continuity of Vv and V a on {x : ||x|| £ 5} 
forces 
(2.41) 


hm F n (t) = Fi(0) = v(x)u ■ (v(x)Vg(x) +g(x)Vv(x)) = u • ( v(x) 2 Xg(x) - g(x) 


VW(x) 


m 


uniformly for (x, u,f) £ {x : ||x|| £ S} x S 2 x [0,T], In conjunction with Lemma [2771 this yields 


sup \Jn ■ E 

(x,u):||x||eS,u£S 2 


F n (x,u,N^ n) {u)) - Fi(x,u,0) (n ( h) ( u)) l {JV W(u)<T} 


< 


sup { F n (x,u,t) - A(x,0) l yjn sup E (x (n) (u)) 

■Ilirll^.Qncsa l ) (x,u):||x||g5,uGS 2 V ' 


(x,u):||x||G5,uGS 

-a o 


as n —> oo. This combined with (12.281) shows 
(2.42) 

hm sup E y/nF n (x,u, N^ n \u)) (N^ n \u)\ l {J v(n )(u)<T} 

n ^°° (x,u):||x||eS,ueS 2 L v 

Combining these calculations with the expansion (12.371) . we have 


g 2 (x)v 2 (x) 


Fi(x, u, 0) 


lim n 3/4 E [l7i7V (n) (U) 


-n 1 / 2 

= lim ———-E 


(2.43) 


n -¥oo 2 g(x)v(x) 
-1 


g(x) 3 v(x) 3 

-1 


E 


U)) 2 1 {jvW(u)<t} 
VW(x) 


17iU • v(x) Vg(x) - g(x)- 


v(x) 2 g Xl (x) - g(x) 


U Xl (x) 


m 


□ 


3g(x) 3 n(x) 3 

This completes the proof. 

Let (X^,7)j")fe>o be a Markov chain with transition operator 

Qnf (Xj t) = E [/ (y [x, w(x)U, N n (x, U )],t + n- 3 / 4 7V"(x, U)) 

as defined in (11.161) . Note that (X%)k>o is a Markov chain with transition operator (11.71) (with 
appropriate scaling), while the second coordinate keeps track of the time between collisions. 
That is, if we start the chain from (x, 0), then n 3 / 4 Tj," = TJf is the time at which the fc’th 
collision occurs. 

Lemma 2.10. Let 

Mx,i) :=nE[(X?,7D- (x,t) | (X£,T 0 n ) = (x,t)] 

be the scaled drift of (X£, Tff)k>o- Then 


lim sup 

n_>0 ° (x,t):||x||eS,t>0 


/z n (x,t) - 


-1 


2VW(x) V.g(x) 


3g(x) 2 [ mv(x) 2 g(x) J ’ g(x)v(x) 


= 0, 


uniformly on {x : ||x|| £ S} x [0, oo). 

Proof. The result for the time coordinate is simply a restatement of (12.271) . By symmetry, it 
suffices to prove convergence for the first spatial coordinate only. By definition of Xf , this is 

nE[X? - Xl \ (X£,T 0 ") = (x,i)] = nE [y n>1 (iV^(U)) - xi . 

Let A s(S) := Aj(,S') as defined in (12.61) . Since A g(S) is increasing in n, we have by construction 
that ||y„(t)|| = r n (t) £ S' for all n £ N and t £ [0, A^S)]. So in particular, y n (t) is bounded 
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nE 

Vn,l (iV^(U)) -Xx 

= nE 

(y n ,i (Y")( U)) - xr) 1 {JV( „)< A , (S)} ' 



+ n E 

(yn, 1 (x (n) (U)^ - *1) l{JV(")>Aa(S)} 


away from dV. To compute jj n .i we first split the expectation on the events {7V( n )(U) < A^S 1 )} 
and {N^( U) > A< 5 (5”)}. This allows us to write 


(2.44) 


By Lemma 12.81 the second term of (12.441) converges to 0 uniformly on {x : ||x|| £ S}. To 
compute the first term of (12.441) . we utilize a second order Taylor expansion of y n ,i(t) evaluated 
at t = (U) which yields 

2/n , 1 (x,AT(")(U))-a ;i =u n (x)f/ 1 iV(")(U) + ^ n , 1 (r) [N^ n \ U) 

for some 0 < r := r(x,7V n (U)) < N n (TJ). Hence 


n E 


(2.45) 


(yn,i (jvW(U)) -* 1 ) l { ^)< Aj(s)} =E nu n (x)t/iY n) (U)l {]V („)< A5(5)} 


-E 

2 


V n ,l(T)N W (V) 2 l {N M<A tm 


Note that 
E 


nt'T.(x)^ • N^(V)1 {nM>As{s)} \ < sup|u(x)|n 3 / 4 E [n^( U)l {Ar( „ )>Ai(s)} 

uniformly for x such that ||x|| £ S as n —> 00 by Corollary 12.61 Thus the limit of the first term 
of (12.451) can be computed by a direct application of Lemma 12.91 


(2.46) 


lim E 

n—> 00 




= lim v(x)E [n 3/4 f7i • iV (n) (U) 


1 ( W Xl (x) g Xl (x.) 


3g(x) 2 \mv(x) 2 g(x) 


uniformly on {x : ||x|| £ S'}. We now compute the limit of the second term on the right hand 
side of (12.45[) . 

Since r(N^ (9)) < N^ n \d), by Lemma T2.31 and Lemma [ 


sup 

(x,u):||x||es,ues 2 


sup 


E 


V^M,i(t(Y">(u))) + 


1 {NM<A s (S)} 


< 


sup 


(x,u):||x||G5,uGS 2 lt€[0,A$(£)] 


y/ny n! i(x,u,t) + 




m 


•E 


\/niV^(u) 2 l{^(„)< Aj (5^ 


as n —> 00 . This together with Corollary 12.61 and Lemma [577] yields 


lim sup 

TWO ° (x,u):||x||e5,u6S 2 


E 


ny n A T (N (n) (u)))iV (n) (u) 2 l {Ar( n)< A5(s)} 


2W Xl (x) 


mg(x) 2 t)(x) 2 


= 0. 


Since this convergence is uniform in u, we can evaluate the limit of the second order term of 
equation (|2.45l) by 


?y nil (r(iV(")(U)))Ar(")(U) 2 i {JV( „ )(u) < Ad(s)} 


If XI ( x ) 


(2.47) lini E / \o / \2 

n ->00 L 2 \ v )- <u mg(x.) z v(x) z 

uniformly on {x : ||x|| G S}. So by adding the right hand sides of l2.46l and 12.471 we have 

lim nE [AT - xi|(X£,7T) = (x,i)] = — 1 ^ l(x) ■ fc(x) 

n—too 

uniformly on {x : ||x|| £ S} x [0,oo). 


3g(x) 2 \ mv(x) 2 C/(x) 


□ 
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Lemma 2.11. Let ((X£, T^))k>o be a Markov chain with transition operator as in (11.1611 and 
let 


Then 


: = nE 


: = nE 


[(XI, - Xi ) (X? tj - Xj) | (Xg , Tq ) = (x,i)] 
[(*"< ^ *,) (71" - t) | (Xg , Tq 1 ) = (x,t)] 
'(7?-t) 2 |(XS,7?) = (x,i)' 


(1) lim 

n—>• oo 

(2) lim 

n—>• oo 

(3) lim 

n—>• oo 


CT "’^ (X,<) 3g(x)2^' 

sup <^,*t(x,t)=0. 

(x,i):||x||e5,i>0 

sup er 2 Ax, t) = 0. 

(x,t):||x||es,t>0 ’ 


sup 

x,£):||x||e5,£>0 


= 0, where <5^ 




Proof. Note that © is immediate consequences of Lemma 12.71 while © follows from ©, © 
and an application of Cauchy-Schwarz, It remains to establish ©• As in the proof of Lemma 
12.101 we have 


(2.48) 




+ nE 


(y n ,i(N {n) ( U)) - a*) {yn,j(N (n) ( U)) 

[(l/n,i(JV (n) (U)) -**) (i/„j(JV w (U)) 


- XjJ 1{NM<A S (S)} 

— X j) ^{NM>A s (S)} 


The second term on the right hand side goes to 0 by Lemma 12.81 To evaluate the limit of the 
first term on the right hand side of (12.481) . it is enough to use a first order Taylor expansion of 
y n(t) 

y n ,i( x > JV (n) (U)) -Xi = y n ,i{T)N^( U) 

for some 0 < r := r(x,7V n (U)) < 7V”(U). Hence 


(2.49) nE [(2 /. m (XW(U)) - ®<) (»„# W (U)) - a*) 1 {sW <a s(S )} 


= nE 


&,iWi,j(T)iV (n) (U) 2 l {JV („)< Aj(s)} 


Note that for t < As (S'), and x such that ||x|| £ S, ||y„(x, u, i)|| £ S'. By the mean value 
theorem and the fact that S7U is bounded on {z : ||z| £ 5"}, there exists C > 0 such that 


sup sup |n 1/4 y n y(x, u, t) - v(x)uj| < 

(x,u):||x||G5,uGS 2 t£[0,A,5(5)] 


mn V 4 z:||z||eS' 


sup U Xi (z)As(S) < 


C 

nl/i 


Since t{N^(9)) < N^ n \d), we can apply Lemma 12.41 and Lemma 12.71 to see that 


lim nE 

n—t oo 


yn,<Wj/n,i(r)JV ( " ) (U) 2 l {ArW < A4(S)} 


lim u(x) 2 E 

n—>■ oo 


UiUj{n 1/2 N^ n \ U) 2 ) 


2 , 
3ff(x) 2 


uniformly on {x : ||x|| £ S}, which completes the proof. 


□ 


3. Convergence of the free path process 

In this section we will prove the convergence of the free path process (X£, 7).") fcgNo with 
transition operator Q n defined in (11.161) . Recall that V 3 = {x £ M 3 : ||x|| £ V}. The following 
is the main theorem of this section. 
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Theorem 3.1. Let (Xt)t >o be a diffusion on T >3 whose generator G acts on functions f £ C 2 (T> 3) 
with compact support in T>\ by 


(3.1) 


G/(x) = -i_A/(x) - 


Vg (x) 2VW(x) 


• V/(x) 


3g 2 (x) 3 5 2 (x) \ g(x) mv(x) 2 

and killed if/when the diffusion TZ := ||T|| hits the boundary ofT>. Consider any l,u £ T>° with 
l < u and start the process (X£,7^") i;eNo at (x, 0), where l < ||x|| < u. Define the stopping times 

r "u := inf{fceN 0 :||X£||^M} 

and 

n,u ■= inf{i > 0 : ||T t || £ [l,u]}. 

Then, as n —> 00 , the family of continuous time processes ^X™ nt j Ar „ , 7j^ t j At « ^ converges 
in distribution on the Skorokhod space to the diffusion 

ds 


Xi 


t/\Ti iU , 


g(X s )v{X s ) 


t> 0 


with initial conditions (x, 0). 

As an immediate corollary we get the following result for the diffusion 1Z. 

Corollary 3.2. Suppose X is the diffusion defined in Theorem \3.1\ Then 1Z := ||«Y|| is a 
diffusion on T>° with a generator G r that acts on functions f £ C 2 (T>) with compact support in 
V° by 

' d r g(r) _ 1 2d r U(r)' 

3 g 2 (r) V g(r) r 


(3.2) 


G,/(r) = 


1 d? r f(r) - 1 


3 g 2 {r) 


2 (r) 


d r f(r). 


Proof. This follows by using Ito’s Lemma to the diffusion X from (13.11) together with the following 
two observations: 

a. If g(pc) only depends on ||x|| = r then 

V.g(x) = d r g{r)~. 

r 

b. If X =: (X\,X 2 , Xf) and (£?i, B 2 , Bf] is a standard three dimensional Brownian motion 
then 

3 

^2 XidBi/lZ 

i =1 

is a standard one dimensional Brownian motion. This follows from the fact that 


Y^XidBi/lZ 


.»=1 


3 X 2 

= V -4 dt = dt. 

^ TZ 2 


i=1 


□ 


We prove Theorem 13.11 at the end of this section. The proof utilizes Theorem IX.4.21 from 
jJS03j which gives sufficient conditions to prove a continuous time step process converges to a 
diffusion. We add this result below in order to make it easier for the reader to follow our proof. 

Theorem 3.3. Suppose that for each n £ N, X n is a pure jump Markov process. That is, its 
generator has the form 

A n f(x) = J[f{x + y)~ f(x)\K n (x, dy) 
where K n is a finite transition kernel on W 1 . Then define b n and c n by 

(3.3) b n {x) = J yK n (x, dy), c n ’ ij (x) = J yYK n (x,dy). 
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Let b, c be continuous functions on R d and suppose X is a diffusion whose generator is given by 
(3.4) Gf(x) = ^2,bi{x)Dif{x) + i c l \x)D i:j f{x) 

i<d l<i }t 7<d 

and defines a martingale problem with a unique solution (see Assumption IX.4-3 from |.TS03| 1. 
Assume that 

(i) b n —> b, c n —> c locally uniformly; 

(ii) sup x . lx \< a f K n (x,dy)\y\ 2 t {M>e} 0 as n f oo for all e > 0; 

(iii) u n —^ v weakly, where v n and v are the initial distributions of Xff and Xq respectively. 
Then the laws L(X n ) converge weakly to P = f P x v(dx), the law of the diffusion process X 
started with the initial distribution u. 


Because we need to stop the diffusion before it enters a domain where the particle’s velocity is 
not bounded away from 0, the easiest way to apply Theorem l3.3l to our situation is to introduce a 
cutoff version of the chains ((XjJ, 7). n ))fc>o, that behaves like this chain on {(x, t) : l < ||x|| < u} 
and obviously satisfies the conditions of Theorem 13.31 outside this domain. To this end, fix 
l*,u* £ V° such that l* < l < u < u* , and let <f> : M —> [0,1] be a smooth function such that <f> 
is identically 0 in a neighborhood of 0 and in a neighborhood of oo, and (f>(r) = 1 if r 6 [l*,u*] 
and let £(x) = (||x|| Vl)Aa. Moreover, define 


b lu (x,t) = b lu (x) 


( - 1 

[2V(#/)(x) V(^)(x)1 

1 ^ 

\3g(C( x )) 2 

_rm;(C(x)) 2 g{ C(x)) 

5(C( X )MC( X ))/ 


and 


c^CM) = c d( x ) = 


3ff(C( x )) : 


;Sijt(i < 3). 


Let A be uniformly distributed on {—1, l} 3 , and define the transition operator R^ by 


i£"/(x,t)=E 


/ 


(x, t) + n _1/2 ^/cft(x)(A, 0) + n-'b 1 



The next two results follow easily from these definitions. 


Proposition 3.4. Let Qi u act on f £ C 2 (R 4 ) with compact support by 

3 lu( \ 4 

(3.5) ft«/(x, t) = Y C? % X Ai/(x, t)+Y &i“( x )A/(x, t). 

1=1 2=1 
The martingale problem for Qi u is well posed. 

Proof. The generator Qi u is slightly degenerate because c l 4A = 0, but this is not a serious compli¬ 
cation, as is easily seen from the SDE perspective. In particular, note that if (Zf, Zf, Z 3 , Zf) t > o 
solves the corresponding SDE 

(3.6) dZ i = y/c l ff(Z\Z 2 ,Z 3 )dB i + b lu (Z\Z 2 , Z 3 ),dt 
for 1 < i < 3 and 

dZ A = b lu (Z\Z 2 ,Z 3 ) 4 dt, 

then Z 4 is a measurable, deterministic function of (Z 1 , Z 2 , Z 3 ) and Z 4 does not appear in the 
SDE for (Z 1 , Z 2 , Z 3 ). Thus we need only establish existence and uniqueness in law for the SDE 
(13.611 for (Z x ,Z 2 ,Z 3 ). However, this is standard since the drift coefficients are bounded and 
continuous and the diffusion matrix is continuous and positive definite, see e.g. (EK861 Theorem 
8.1.7], ' □ 

Proposition 3.5. Let (Z%)k>o be a Markov chain on R 4 with transition operator R„. Then, 
uniformly on R 4 , we have 

(1) lim nE(Z” - z | Z£ = z) = b lu { z) 

n—foo 

(2) lim nE((Z 4i - Zi )(Z^ - Zj ) \ Z£ = z) = c&(z) 
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(3) lim max||Z5J +1 — ZjJ|| = 0 almost surely. 

n—>oo k >0 

To handle the cutoff, we will need the following elementary result about convergence of hitting 
times with respect to the Skorokhod topology. 


Lemma 3.6. For 7 G D[ 0, 00 ) define r x ( 7 ) = inf{£ : 7 (t) > a;} to be the first passage time of 7 
across x. Consider an f € C\ 0, 00 ) and an a > 0 such that 

0 < T a (f) < 00 and inf {£> r a (/):/(£)> a} = t 0 (/). 

If(fn) n >i is a sequence in D[ 0, 00 ) such that f n f in D[ 0, oo), then T a (f n ) — ► r a (/). 


Proof of Theorem \S.1\ Let us introduce the cutoff transition operator defined by 


7) J to,) 

1 Rn(^t) if ||x|| i [l,u] 


Let ("x[“’ ra , 7 ^“’") be a Markov chain with transition operator Q^/(x, t) started from (x, 0 ) 
V / fc>0 

with ||x|| e (l,u) and let r/n = inf{fc : ||X)“’ n || ^ [Z,u]}. By the definition of the cutoff transition 
operator, we have the identity in distribution 



so it suffices to prove the result for the cutoff chain T^ U 'J . In order to apply Theorem 

13.31 directly, we Poissonize the chain. That is, let (Tt)t>o be a rate one Poisson process and 
consider 


fvlu.n '~rlu,n (j.\\ _ (-yj-lu,n ^r'lu,n\ 

[t)j WJ t > 0 — [^r t > 'r t ) t>0 • 

Since T nt /n —> t uniformly on compact sets almost surely, it is enough to prove the result for 
(X. lu,n (t),T lu,n (t)) t>Q , where the relevant stopping time is now 


rf n = inf{i: ||X'“'"(t)|| * Ml- 

Note that (X iu ’ n (n£), T lu ’ n {nt)) t>Q is a pure jump Markov process with generator 


g l *f{x,t)=n(Q l *-I)f{*,t), 

see e.g. |Kal021 Proposition 19.2]. For this process, condition (i) of Theorem 13.31 is satisfied with 
b = b lu and c = c lu by a direct combination of Proposition 13.51 Lemma 12.101 and Lemma 12.111 
That condition (ii) of Theorem 13.31 is satisfied is an immediate consequence of Proposition 13.51 
Lemma 12. 11 1 ( for the time component), and Lemma 12.81 where to apply the last Lemma we note 
that, because velocity is bounded away from 0 and 00 on [l,u], the condition that N^ > e in 
Lemma 12.81 is comparable to the condition that the size of the jump be at least e in condition 
(ii) of Theorem 13.31 Consequently we have 

(X'“'"(nt) ) r , "-"(n*)) t > 0 =► (Tf, 7?% 0 


where (A^ u , T t lu ) t>0 is a diffusion with generator (13.51) started from (x, 0). Using the Sko¬ 
rokhod representation, we may assume this convergence happens almost surely. It follows from 
Girsanov’s theorem that ||<T iu || satisfies the hypotheses of Lemma HTTijl and, as a result, this con¬ 
vergence happens jointly with the convergence of rff/n to rf u = inf{£ : ||T’ Z “|| ^ [/, u]}. Therefore 

(X‘“'»((ni) A rff), T lu ’ n ((nt) A »£))*>„ =*► (<V“^AV“) i>0 

as n —> 00 . From the SDE perspective, see Proposition 13.41 it is clear that the generator for 
((T) z “, T t lu ))t>o agrees with the generator for 


(' XU 1 9(WX S ) 


t> 0 
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on smooth functions /(x, f) with compact support in {x : ||x|| £ [Z*,tt*]} x R, and thus we have 

ds 


( \slu ‘-T'lu \ _ IV 

™'tf\r\ lu ’> ' t/\r} lu ) d l 


tAT! '“’7 0 g(x s )v(x s ) 

As a result of de-Poissonizing the process and (13.81) , we have 

ds 


t> o 




g(X s )v(X s ) 


□ 


t>o 


4. The Process on its natural time scale 


In this section, we study the convergence of the full trajectory of the particle. To reconstruct 
the full path from the free path process, we also need to keep track of the direction of reflection. 
To this end, we need to keep track of the reflection times and the direction of reflection. That 
is, we will look at the Markov process (XJJ, 7^", U^) fcgNo started from (x, 0, 0) with ||x|| £ V and 
transition operator 


Qnf (x, f, /l) = E 


f (y [x, t(x)U, IV" (x, U )},t + n- 3 / 4 iV"(x, U), u) 


By (11.41) the full trajectory ( X. n (t)) t>0 of the particle with energy E n moving in potential U n 
and scattering density g n is then constructed by 


(4.1) 


x"(t) := y n (Xfc, u(X£ )UJJ , t — n 3 ^) 


for t £ [n 3 / 4 T fe n ,n 3/4 7; n +1 ). 

The following theorem is the main result of this section. 


Theorem 4.1. Let (X"(f)) denote the full trajectory of the particle as defined in (14.11) and let 

% u ■= inf{t : ||X"(t)|| ^ 

For any fixed l,u £ T>° with l < u, as n —> oo we have the following convergence in distribution 
on £)(R + ,R 3 ): 

(x”((n 3 / 4 f) A ^ u )) t>0 (X(m A r ; , u )) t > 0 
where X is as in Theorem no and fl is the time change given by 

m ■■=!([ (t) 

\J o A(T(s A Ti tU j)) 

and I is the inverse operator defined by X(/)(f) = inf{s : /(s) > t}. 


To prove a limit theorem for the full path, we must invert the time process {T^ l )k>o- This is 
easiest to do if we extend this process to [0, oo) by linear interpolation: 

(4-2) T n (s) := TJj + (s- LsJXTJj +1 - TJj), Vs £ [0, oo). 

Defining X (f)(t) = inf{s : f(s) > t} we have 

(4-3) X”(n 3//4 7^”) = X^ )(7?)J = X£. 

Consequently, we start by examining the convergence of ^X" I ( 7 - n )( t )j ^ , and then prove that 

filling in the true values of the path on the intervals does not affect the limit. 

The cutoffs do affect the time change in a technical way. In particular, the stopped process 
(' T n (t A Tf u ))t> o does not have invertible paths. To accommodate this, we use the time change 
fl n : R+ —> R+ defined by 


(4-4) 


0„(f) = X 



+ (•- 


1 


A(X"(tA )) 


(t), 
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where A(x) = g(x)u(x). The choice of the particular linear drift is because Theorem l3.ll together 
with the result in the proof about convergence of hitting times, implies that asn->oo 
(4.5) 

( X Lntj AT-y^AT^ +(t~ n l7 -”u) + A (xn( T n ))) =* / A(* sAT ,j) 


t> 0 


t>0 


The next result allows us to invert the time process. 


Lemma 4.2. If f £ .D(R+,R+) is continuous and strictly increasing with lim^oo f(t) = oo, 
then 1(f) £ D(R+,R+) and I is continuous at f. 

Proof of Theorem \f-l\ By (14.51) and Lemma 14.21 we have 


(4.6) 


Ar" u ) ^™(^) t>0 ^ (^tAri.u) ^(t)) 


t >0 ' 


Note that the process (fl(f) t >o is supported on C(R+,R+) and, consequently the composition 
map (/, g) H > f ° g from Z?(R + ,R 3 ) x C(R + ,R_|_) —> D(R + ,R 3 ) is continuous, see e.g. lBi!991 p. 
151]. This combined with (14.61) implies that 


(4.7) 


(’ 


dX(T“)(t)J At 


) = ( X 

/ t>o V 


|n£2 n (t)J At 


) t>0 ^ (^(t)An,J 


t> 0 


as n oo. For T > 0 fixed, by SSD , we have that 


sup 
o <t<T 


Note that 


(n 3/4 - x y (T ., (l)J 

< sup 

k<\nn n (r/\Tj^ 


sup 

T£_ x <n-VH<T£ 




P ( sup 

\0 <t<T 


< e + P sup sup 

\k<n 2 T£_ x <n~ 3 IH<T£ 


> e 


n^n (T A 77AJ < nSl n = t£ u = 0(n) 

in probability since n _1 T™ u —> n iU < oo a.s. Thus if we fix e > 0, for n sufficiently large we have 

X" (n 3/4 (t A 75 J) - Xy (T . )(t)J AT . 

y„ (x£_ 1; U k - u t- - Xj^ t {k < T n +1} > e) . 

TT V 7 / 

The probability on the right hand side can be made arbitrarily small by observing that by, 
Corollary O for C > 0, 

P(n 3/ \T k n - Tff_ x ) > C , |||X£_ 1 || S [l,u]) < sup P(W<">(x,u) > C) = o(n k ) 

(x,u):||x||e[i,u].u6S 2 

in probability for every k > 1 and by Lemma 12.31 if n is large enough then 

sup ||y n (x, u, t) — x|| < e. 

(x,u,£)e{x:||x||G[i,ti]}xS 2 x [0,C] 


Therefore 


00 : 


( X ” (» a/4 (f A ^".))),> 0 ^ 

as n —> 00 and the fact that 


as n —> 00 follow from path continuity and observing that if u < . 


□ 
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5. Classifying the boundaries of V 


In this section we give conditions under which the boundary points of T> are inaccessible, 
meaning the boundary points cannot be reached in finite time. These results can then be used 
to remove the stopping at u and/or l in Theorem 11.41 

Suppose we have a regular diffusion X with state space the interval (£, r). Every diffusion has 
two basic characteristics: the speed measure m(dx) and the scale function s(x). 

We assume the infinitesimal generator G : T>(G) i-a C&(/) of X is a second order differential 
operator 

Gf(x) = ^ 2 ( x )9 xx f(x) + p(x)d x f(x) 

where a, g, £ C(I) and er 2 (a:) > 0 for all x £ I. Let B(x) := f* 2a~ 2 (y)g(y) dy for some arbitrary 
(fixed) a € I. Then it is well-known that 

• The speed measure is absolutely continuous with respect to Lebesgue measure and has 
density 

m'( x) = 2cr _2 (x)e B ^ 

• The scale function has density 

s\x)=e- B G). 

The domain T>(G) consists of all functions in C&(/) such that Gf £ Cb(J) together with the 
appropriate boundary conditions. 

The boundary point £ is called accessible when 


and inaccessible when 


w; m'(rf) dr/^j s'(y) dy < oo 

/ (/ s '^ d y = °°- 


Similarly, one can classify the boundary r. 

As we have shown above in Corrolarv l3.2l the one-dimensional diffusion defining the limiting 
radial process has generator Q r that acts on compactly supported functions in C 2 (V°) as 


(5.1) 


Qrf{p) ■= Pr{p)f\p) 


1 {P) 


f"{p) 


g'{p) , 1 d r u{p) 


3g 2 (p) \ g(p) p (E-u{p)) 


f'(p ) 


3 g 2 {p) 


f"(p )• 


If a £ V° then the density of the scale function will be 

rV 2 Hr(p) 


( r y 2p r (p) - 

s(y) = exp l“A ^M dp . 


= exp 


r v (g'{p) 1 d r U{p) 

a \g(p) P ( E-U{p)) 


dp 


= exp (In (g(y)) - ln(y) - In (E -U(y)) - C) 

g{y) 


= exp In 


y{E~U{y)) 


-c 


a{E-U(a)) 

g(a) 


g(y) 


d y(E-U(y)), 

where C := In (^ a (^E 9 -u(a)) ) ■ a result, if we fix an arbitrary c £ V° the scale function will be 
given by 


(5.2) 


s{x) = / s'{y)dy 

J C 


a{E-U{a)) 

g(a) 


i: 


g(y ) 


y(E-U(y)) 


dy. 
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The speed measure density for the radial diffusion will be 

2 2 g(a) 


(5.3) 


m'{x) = 


<jf(x)s'(x) a(E — U(a)) 
One can then find the speed measure by setting 

2 g(a) I 

a(E — U{a)) J, 


i(J) := J m'(x)dx 


g(x)x(E -U{x)). 


g(x)x(E —U(x))dx 


for any Lebesque measurable J C D°. 


Remark 5.1. If s(£) = —oo then I cannot be reached in finite time and is therefore inaccessible. 
Similarly, if s(r) = oo then r is inaccessible. 

Proposition 5.2. Let D° = (h_,h + ). 

i) Suppose that there exists a > 1 such that when y \.h- 

1 


= 0{(y-h-)~ a ). 


y(E-U(y)) 

Then h— is inaccessible. 

ii) Suppose that there exists a < 1 such that when y J. h- 

1 


= 0((y-h-)~ a ). 


y(E-U(y )) 

Then h- is accessible. 

In particular, if h_ =0 and |W(0)| < oo then 0 is inaccessible. Analogous results hold for h + . 


Proof. The main parts of the proposition follow easily from (15.211 and Remark 15. 11 
If h- = 0 and |W(0)| < oo then, since by assumption E > U(0), we have 

1 


y{E-U{y)) 




□ 


5.1. Constant force. Suppose we have a constant force field 

U(x) = Cx 

and V = [0, h] for some h > 0. We know that the speed has to be zero at h so 

E — U{h) = E — Ch = 0. 

This forces E = Ch and thus 


1 


1 


y{E-U(y )) yC(h — y )' 

By Proposition 15.21 both 0 and h are inaccessible boundary points. 

5.2. Newtonian gravity. We consider the case when the potential is that of a gravitational 
force directed towards the origin. That is, the potential function is given by 


U(p) = — 

P 

for some constant k > 0. Suppose the the total energy of the particle is positive E > 0. We have 


s(z) = 


(.Ea + k)a 

9(a) 


g(y ) 


(Ey + k)y 


dy. 


Since g is bounded above and bounded away from 0 on R + , for y J. 0, 


g(y) 


\/(Ey + k)y 


= o\- 
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which implies lim^o s(x) = — oo. For y t oo one has 


so that lim^oo s(x) 


g(y) =o(-) 

(Ey + k)y \y 2 ) 


oo. By Proposition 15.21 the boundary points 0 and oo are inaccessible. 


6. Higher dimensions and general potentials 


Although our model was motivated by the three dimensional Lorentz gas, the process with 
uniform reflection direction makes sense in any dimension and our methods apply essentially 
without change. In this section we explain how our results could be generalized to higher 
dimensions and to non spherically symmetric potentials. Consider a particle in for rf £ N 
under the influence of a force field with potential energy U{x). Let S’ d_1 be the unit sphere 
in R rf . In this model upon reflection at x = (xi, ..., Xd), the particle starts its path in the 
uniform direction u = (ui,... ,Ud ) £ S' d_1 and then travels under the influence of U(x). The 
only difference is that some of the constants change because we are integrating u over S d_1 
instead of S 2 . Using the fact that if u = (ui ,..., Ud) is uniformly distributed on S d_1 then 
Em 2 = h f° r all i £ {1,..., d} yields the following results. 

Let ((X£, T)?))k>o be the free path Markov chain, where is the location of the particle at 
the time of the fc’th reflection and is the time of the fc’th collision. 

Theorem 6.1. Let 


■■= nE[(X?,77*) - (x,i) | (Xg, Tq ) = (x,i)] 

be the scaled drift of (fX. 7 ^,Tf l )k>o, a-nd let S be a compact subset ofT>°. Suppose the equivalents 
of Assumptions ( A 1 )-( A 5 ) are satisfied. Then 


lim sup 

n ^°° (x,t):||x||eS,i>0 


/x n (x,t) - 


-1 


(d-l)VW(x) Vg(x) 


<?(x) 


1 


dg(x.) 2 mvfx) 2 
uniformly on {x : ||x|| £ S} x [0,oo). Furthermore, let 

< i3 -(x,i) := nE [(XI, - Xi ) (X” . - Xj ) | (X£, T?) = (x,i)] 


g(x)u(x) 


= 0, 


< it (x,t) := nE [(XI, - Xi ) (7? - t) | (X 0 ",7^) = (x,i)] 


<t(x,t) :=nE (T? -tf | (Xq , Tq ) = (x,i) 


Then 


(1) 

lim 


sup 

2 

°n,ij 

(x,i) - 


n—foo 

(x,t) 

l|x||€5,i>0 


(2) 

lim 


sup cr 2 ( 

x,t) = 


n—foo 

(x,t) 

l|x||G5,t>0 



(3) 

lim 


sup a 2 t (x,t) = 


n—foo 

(x,t) 

||x||G5,£>0 




d- 1 


dg (3 


\2 W V 


= 0, where 5,j = t{i = j }. 


Theorem 6.2. Let (AF(t),t > 0) be a diffusion on T> C whose generator 5? acts on functions 
f £ C 2 (V) with compact support in C 2 (V°) by 


&f(x) 


v(x) 

dg(x) 


A/(x) 


v(x) ( jd-l)XU{x) Xg(x) \ 
g(x) \ dmv 2 {x) dg{x) ) 


Let (X”(f),f > 0) be the trajectory of the particle and assume the equivalents of Assumptions 
( A 1 )-( A 5 ) above hold. Fix a set J C T>° and define 


t n j = inf{f > 0 : X"(t) £ J} 

and 

tj = inf{t > 0 : 3£(t) ^ J}. 
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Then as n —> oo we have the following convergence in distribution 

(x"(n 3 / 4 tAi}),t>0) -A (&(tATj),t>0). 

Futhermore, if the boundary of T> is inaccessible then we can remove the stopping in the above. 

Acknowledgments. We thank Simon Griffiths, Cristina Costantini and Steve Evans for 
very helpful discussions. 
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